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Abstract 

We have discovered conjectural near-addition formulas for Somos sequences. We have preliminary evidence 
suggesting the existence of modular theta functions. 

1 Introduction 

In 2001, as part of a research project investigating more efficient public key cryptography (PKC), Rich Schroeppel 
asked Bill Gosper to look for Somos sequence addition formulas. Gosper found some very interesting results 
immediately, and further developments continued through 2003. Cheryl Beaver and Schroeppel also investigated 
modular versions of the dilogarithm function Li2(x) and the Trilogarithm Lis(o:) [6], and Schroeppel did some 
preliminary work on modular theta functions. 

Somos sequences and theta functions are both promising approaches for use in public key cryptography. 
Cryptographic applications of Somos sequences are explored in pQ. Our dilogarithm and trilogarithm results are 
interesting, but it's not obvious how to apply them to PKC problems. 

Further development of modular versions of the special functions of numerical analysis seems possible. Likely 
candidates are the error function, the logarithmic integral, the gamma function, and perhaps Bessel functions and 
hypergeometric functions. 

The main stumbling block is that inequalities, limiting processes, and infinite series are unavailable, and we 
must fall back on functional equations for much of the work. Formal differentiation sometimes works. Functional 
equations are very limited for the error function, but more variety is available for the other special functions. 
We have not explored the p-adic possibilities, which might permit the reintroduction of some of the forbidden 
concepts. 

Both Somos sequences and theta functions have near-addition formulas: equations that relate f(x + y)f(x — y) 
to f{x) and f(y) and / of nearby x and y values. These can be used with the well-known double-and-add method 
to calculate function values at large multiples of the argument. 

Gosper's results on near-addition formulas for Somos sequences are reported in section 2. Schroeppel has been 
able to prove a few of the formulas, including two of the determinant identities for Somos4 Q]. Section 3 details 
our brief excursion into modular theta functions. There is no conclusion: this research seems open ended. 



*rwg@osots.com 

trschroe@sandia.gov Sandia is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company, for the 
United States Department of Energy under Contract DE-AC04-94AL85000. 



1 



1.1 Somos Sequence Background 



RS first learned of Somos sequences from Michael Somos, around 1988. The Somos sequence of order N begins 
with a block of N Is, and is generated from a simple non-linear recurrence. For Somos4, the recurrence is 
akdk+4 = a.k+iak+3 + a k+2 ; for Somos5, it is dkCik+5 = ak+iak+4 + dk+2Cik+3', Somos6 and 7 follow the patterns for 
4 and 5, with more terms in the folded dot product. The recurrence can also be used to extend the sequences in the 
negative direction; they are palindromic. Somos4 begins 1,1,1,1,2,... . One surprising property is that all the terms 
in Somos4 - Somos7 are integers. This was discovered back in the 1940s by Morgan Ward 9, 10 for a sequence 
including Somos4. He called his sequences Elliptic Divisibility Sequences; Somos4 is the odd numbered terms from 
a particular EDS. RS was at an MSRI number theory workshop shortly after learning of the sequences, and the 
group spent some time trying to prove integrality. Eventually Dean Hickerson and Janice Malouf independently 
proved that Somos6 is integral. Our experimenting showed that we could modify the sequence initial values in 
various ways while apparently keeping the integrality property. Also, introducing algebraic coefficients into the 
equation, such as ak(ik+4 = xak+ia k +3 + ya\ +2 , often produced polynomials with integer coefficients, rather than 
the expected ensemble of rational functions. (Of course, this doesn't matter much when the values are interpreted 
mod P, or in a finite field, which can handle fractions just fine.) The raw integer values of the sequences seem to 
grow roughly as C . The Somos4 and Somos5 sequences have a close connection with elliptic curves and classical 
theta functions. The higher order Somos sequences may be connected to hyperelliptic curves. There's a moderate 
amount of background material scattered around the net; Jim Propp's Somos page [7], and Sloane's sequence 
database [S], and Zagier's problem 5 Q2], have useful material. Background on theta functions is available in the 
Abramowitz & Stegun [2] compendium of special functions, now available on the net (as a scanned photocopy) at 
®. 

2 Somos Sequence Near- Addition Formulae 

Summary: Somos(n:r) is calculable in O(logn) time from three values near Somos(:r), at least for orders 4 and 
5. Orders 6 and 7 require longer intervals of values. Along the way, we find addition formula? for Somos and 
Somos-like sequences of polynomials and algebraics, and reduce some fifth order recurrences to fourth and third 
order. We find three-term, four-variable relations for most of these, as well as for ordinary i? functions. A sequence 
of polynomials obeying the Somos4 recurrence has a particularly nice doubling formula. Many of these results 
fall out of a very general determinant identity. For certain algebraic "Somos" sequences, we find closed forms in 
terms of Chebychev polynomials. 

Definitions: a n := Somos4, b n := Somos5, . . . ,e n := Somos8, i.e., 

, 2 

an := ° W ~ 1 "~ 3 "~ 2 = »3-» = ...,2,1,1,1,1,2,3,7,23,59,314,1529,... 

b n ~ b n-lbn-4 + b n „ 2 b n - 3 = = ^ x % ^ ^ ^ ^ _ _ 

bn-5 

Cn — lCn — 5 + C n — 2C n — 4 + c n — 3 oiiiiiiocnno^r 

Cn := = c 5 -n = ...,3,1,1,1,1,1,1,3,5,9,23,75,421,... 

Cn-6 

d n -\d n -e, + d n -2d n -5 + dn-sdn-4 oiillliioEnl'7/1110'7 

d n := = d 6 -n = ...,3,1,1,1,1,1,1,1,3,5,9,17,41,137,... 

dn-7 



(all appropriately palindromic) where the tabulated values start with subscript n = — 1 to show the center of 



2 



symmetry. They are integer sequences until 
420514 

Somos8: ei7 = — , so we're not too interested in Somos8. On the other hand, Somos6 satisfies 

— C„_iC n _8 — C n _2C n -7 + C„_3C n _6 + 34c n -4C n _5 



C n ~9 

which is pretty much a Somos9. And, for all t and it, Somos4 satisfies 

_ (t — 7) a.n-ian-7 + (u — 5t + 31) a n -2a n -% + (4t — u + 1) a n -3a n -5 — ua n -4 

a n -8 

a double continuum of quasi- Somos8s. 

Furthermore, for all t, the sequence s„ := a n satisfies 

_ (6 - t) s„- lSn - 7 + (5t - 130) s n -2S n - 6 + (749 - 4t) s n - 3 s n -5 + (20t - 4) s n _ 4 

Sn-8 

Change of variable: A Somos sequence may be multiplied by any constant. A Somos sequence multiplied by 
an arbitrary geometric progression satisfies the same recurrence, but usually loses its palindrome symmetry. The 
"odd" (Somos5 and Somos7) sequences may also be termwise multiplied by any number of factors of the form 
tan(x + rwr/2) without even disturbing the palindrome property. 

There is, however, no sharp dichotomy between odd and even, since Somos4 satisfies the quasiSomos5 (Quasi- 
modo?) (odd) recurrence 

5an-3ln-2 — a n _4Ct n _l 



as does a n tan(x + ivn/2), etc. 

2 

The sequence s n := r n a n satisfies 

2 

while s n := r n b n satisfies 
Similarly, s n := r" 2 c n satisfies 



O.n-5 



_ r 6 s n - 1 s n - 3 + r 8 s n „ 2 

Sn — ; 
■Sn — 4 



_ r 8 S n -lS„-4 + r 12 S n -2S n -3 
Sn — - 



Sn — 5 



V Sn-lSn-5+r 6 S n -2Sn-4+r 1 S n _ 3 



Sn — 6 

2 

while s n := r n d n satisfies 



_ r 12 Sn-lSn-6 + T 2 ° S n -2Sn-5 + T 24 S n -'jSn-4, 
Sn — ■ 



Sn-T 

Notation: 



_ / Xi, X2, ■ ■ ■ , x n \ , r -I 

D s : = det s x „ +s 

\Vl, 2/2, ■ ■ ■ , Vn ) 
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s xi—y 1 Sx 1 +y 1 S xl —y 2 S xl + y2 
s x 2 -yi s x 2 +yi s x 2 -y2 S x 2 +y2 

s xn~yi s x-n+yi s xn—y2 s x n +y2 



s xi-y n s x 1 +y„ 
S X2 — yn S X2+yn 



S Xn — yn S Xn+yn 



Note that each term of the expanded determinant will have subscripts summing to 2x\ + 2x2 + ■ ■ ■ + 2x n . This is 
decidedly not symmetrical in x and y, so that an identity involving a D operator may yield a new identity under 
interchange of the x and y vectors. 

Conjecture 4: the determinant 



D a 



x, y, 



a u — x &u-\-x Ctu — y Q<u + y Ctu — z a u +z 
0/v — x &v-\-x &v — y &v-\-y &v — z ^v-\-z 
®w — X Qw-\-x — y Qw^-y &w — z 0/w-{-z 



where u,v,w,x,y, and z are arbitrary integers. E.g., 



= 0, 



n-2, 
0, 



0, 1 

1, 2 



-3 On 

2 
1 



-4 dn 

3 
2 



the defining recurrence for Somos4. 

Note that the determinant also vanishes for at := sint, for arbitrary complex u, v, w, x, y, and z. More interestingly, 
experimental Taylor expansion at q — plus several numerical experiments suggest that the same goes for 
at := fij(t,q). The published addition formulae mixing two or more different j are merely the result of choosing 
v,w,y, and z to be things like it/2 and ttt/2 (and 0). (Whittaker & Watson, crediting Jacobi, list numerous 
special cases, suggesting that the more general formula was not yet known.) 

Still more generally, 



= 



fi s (x — u,q) fi t (x + u,q) fi s (y- 
■&s (x - v, q) fit (x + v, q) fi s (y ■ 
fi s {x-w,q)fi t (x + w,q) fi s (y - 



u,q) fi t (y + u,q) fi s (z ■ 
-v,q) fi t (y + v,q) fi s {z 
w,q) fit (y + w,q) fi s (z - 



- u,q) fit (z + u,q) 
-v,q) fit (z + v,q) 
w,q) fit (z + w,q) 



E.g., putting s = \,u = x,v = y 
= 



■&1 (x - w, q) fit (x + w, q) fi-i (z - y, q) fi t (z + y, q) 
- fix (y - w, q) fit (y + w, q) fix (z - x, q) fit (z + x, q) 
+ fii(y- x, q) fit (y + x, q) fix (z - w, q) fit (z + w,q), 



(4vars) 



a three term identity in four variables. 
Conjecture 4.5: The determinant 



D s 



it + 1/2, w + 1/2, iw + 1/2 
a; -1/2, y-l/2, z - 1/2 



Sl+u — x $u-\-x 
Sl+u — x $v-\-x 
S\-\-w—x Sw-\-x 



Sl+it — y Su-\-y 
Sl-\-v — y $v-\-y 
£l-\-w — y Sw-\-y 



where and u,v,w,x,y, and z are arbitrary integers. E.g., 

' n ■ 



= D b 



5 

2 ' 



1 

2 ' 



1 

2 ' 
3 

2 ' 



Sl + U — Z Sll-{-Z 

Sl+u — z Sv+z 
S\-\-w—z Sw+z 



= 0, 
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bn-3 b n -2 b„-4 b n -\ b n s b„ 
12 3 
1 1 2 



= — b n -5 b n + b n -4 &n-l + &n-3 &n-2, 



the defining recurrence for Somos5. 

Somos4 addition formula?: (See the section "Somos4oid polynomials" for a sequence s n with much nicer 
addition formulae than those derived here for a n .) 

Suppose we have four consecutive values a x -i, a x , a x +i, a x +2- Choose s := a, u = x, y = 0, z = —l,v = Q,w = 1 
to get 



D a 



x + 1, 0, 1 
x-1, 0, 1 



D a 



x + 1/2, 1/2, 3/2 
x-1/2, 1/2, -3/2 



fflifl2i a x a x +i a x -ia x +2 
ai- x a x ao ai a-ia2 
0*2-x o, x +i a l a 2 0,0 a 3 

a.2 X a x a x +i a x -ia x +2 
a x a x +2 1 2 
al+i 1 1 



giving us 02x- Alternatively, 

'a; +1/2, 1/2, 3/2 



Da 



2-3/2, 3/2, -1/2 



d2-x£lx-l (1-1(12 «0 <Jl 

ffl3-a; Ox ^0 013 ffll (22 

Oi+i 2 1 



= 0, 



giving us a 2 x-i- 

Now run the Somos4 recurrence one step forward to get a x+ 3 and replace x by x + 1 in the preceding two 
determinants to get the four consecutive values a,2 X -i, &2 X , 02^+1, a,2 X +2- So we can double a nx to «2 ra . 

Now suppose that we have the four values around a nx and also around a(„+i) x . Then 
D a 



nx + 1, 0, 1 
(n+ 1)2-1, 0, 1 



, (n + l)a; + 1/2, 1/2, 3/2 
" ' m- 1/2, -3/2, -1/2 



Ox+lO(2n+l)x »(»+l)i-l"(ii+l)i+2 0(n+l) x 1(„+1) a,+l 

flnill-ni 1-1(12 ClOll 

12-nifln i+l 10 0,3 ai 0,2 

«i+l"(2n+l)i 1(n+l)x-l 0,(n+l)x+2 Cl(n+1) x 1(n+l) x+1 

x x + 2 2 1 

a n X + l 1 1 



= 0. 



So from d„ and d(„ + i) x we get a2„ x and a(2 n +i) x - Thus we can multiply by maintaining eight values. E.g., 
105x <- (532, 522) <- (272, 262) <- (142, 13a;) <- (72, 62) <- (42, 32) <- (22, 2). 
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In principle, we need only maintain two sets of three values, a n x-i, a nx , a„x+i, and 
a(n+i)x-i, a( n +i)x, a,(n+i)x+i, by virtue of the third order relation 

Conjecture 4a ("derived" below): 

2 2 3 3 2 2; 

a x-i a x+2 + a x a,x+2 + CLx-i a x +i + a x a x+i = 4 cix-i a x a x +i a x +2, 
with which we can eliminate a x +2 from: 



<l2x-l Obx-l Clx+2 dx CLx+l 

dx-1 Clx+1 2 1 



d2x-l ~ CLx-1 (li-l ~ «i) Clx + 2 + dx li + l (li-lfli+1 ~ 2ttj] 



a% 1 1 

to get 

2 / 2 2 2 4 \ 

«i-ii2i-i-«i (2a I ,_ 1 a x+1 — cix-i a x a x +i + a x j a,2 X -i 

( 4 4 ,3 23 , o 2 42 c 6 i n 8\ n 

+cia:+i [a x -i a x +i — 4a a ._ 1 a x a x+1 + oa x _ 1 a x a x+1 — t> a x -\ a x a x +i + 2a x j = U. 

Similarly, 

32 / 2 \ / 2 2 2 4 \ 

"i-i^i + a i (2 a x -i a x +\ — a x j [a x -i a x+ i — 5 a x -\ a x a x +\ + a x j ci2x 

,3 / 4 4 3 23 i or, 2 42 6 , o 8\ n 

+Q-X+1 {a x -i o-x+i — oa x _ 1 a x a x+ i + A\Ja x _ 1 a x a x+1 — 14a x -i a x a x+ i + 6a x j — U 

and 

5 2 /o 4 4 1 o 3 23 l00 2 42 n 6 . 8\ 

a x-i a 2 X +i — a x [aa x _ 1 a x+1 — L8a x _ 1 a x a x+1 + ll a x _ 1 a x a x+1 — y a x -\ a x a x +i + a x ) ci2 X +i 

i 5 (a 4 4 io3 23 on 2 42 n 6 i t 8\ n 

+0^+1 [4a x _ 1 a x+1 — 12a x _ 1 a x a x+1 + 2Ua x ^ 1 a x a x+1 — lbo x _i a x a x+l + 7a x ) = 0. 



The square roots plus the size of these expressions probably render them "too cumbrous to be of any importance,' 
but the even coefficients may pay off in some finite fields. 

"Derivation" of Conjecture 4a: By Conjecture 4.5, 



fe/2-3/2, fe/2-1/2, fe/2 + 1/2 
fc/2-1/2, fe/2 + 1/2, fc/2 + 3/2 



fl-iat-2 a-2dfc-i 0-3(1* 

OO C-lOi 1-2 %1 



2 ((4a fe _ 2 a*: — 3afc_i) a fe+2 — 3a fc _ 2 cifc+i + 8a fc _i a k a k +i — 7a|). 



Da 



fe/2-1/2, fe/2 + 1/2, fc/2 + 3/2 
fc/2-3/2, fe/2-1/2, fe/2 + 1/2 



03 dfc (22 dfc+l Ol afc + 2 
2 



(a fc _ 2 ifc — afc-i) cife+2 — flfe-2 flfe+i + 3 a fc _i a fc a fc+ i — 2a|. 
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Eliminating a k -2, 

(22. 3 3 2 2 \ 

a>k-i a k+2 — 4a k -i a k a k+ i a k+ 2 + a k a-k+2 + a k -i a k+1 + a k a k+1 J = 0, 

as desired. 

Nonstandard initialization: You might wonder how this third order recurrence can compute a fourth order 
recurrence with four initial conditions (ao, . . . a-j = 1). First of all, given the palindrome condition and scaling, 
there is only one degree of freedom. I.e., in general, we have 

o_i,oo, • • • = p 2 +p, 1,P,P, 1,P 2 +P,P 2 + p+-,p' i + 2p 2 + 2p+-^- + l,... . (lppl) 

p p 2 - 

When p is a root of unity, the denominators remain bounded and can be scaled out, e.g., 

. . . ,i - 1, l,i - 1,-1, i - 2, 2 - 3i, -1, 13 i + 3, -16 i - 15, -19 i - 44, . . . 

. . . , i (\/2 + l) + 1, \/2, i + 1, i + 1, \/2, i (y/2 + l) + 1, i V2 + 2, i (y/2 + 3) + s/2 + 1, 
i (3V2 + 7) - V2 - 3, i (5^2 + 12) + 4V2 + 2, i (l5-\/2 + 13) - 22V2 - 31, i (8V2 + 8) - 43\/2 - 76, 

i (34^2 + 57) - 190^2 - 287, . . . 

. . . , 2 i \/3, 2, i V3 + 1, i V3 + 1, 2, 2 i V3, i V3 + 1, 3 i V3 - 1, 
-10, 6i \/3 - 8, -21 i V3 - 9, 35 - 9 i V3, 136 - 66 i V3, . . . 

...j(V6-V2 + 2) + (V2 + l)V6 + V2,4,i(V6-Vty 

i (VE - V2 + 2) + (V2 + 1) VE + V2, 2 i + (V2 + 2) VE + 2V2, 

i (2VE + 2) + (3V2 + 2) v / 6 + 4\/2 + 4,i ((4\/2 + 3) \/6 + 11 \/2 + 6) + (5\/2 + 5) \/6 + 15\/2 + 20, 
i ((5^2 + 6) \/6 + 6\/2 + 10) + (15^2 + 18) V6 + 38\/2 + 46, 
i ((46\/2 + 67) VE + 107\/2 + 156) + (54\/2 + 89) VE + 131\/2 + 212, 
z ((21(h/2 + 31l) \/6 + 523\/2 + 772) + (250\/2 + 34l) \/6 + 635\/2 + 860, 
i ( (963\/2 + 1410) VE + 2346\/2 + 3434) + (l383\/2 + 1934) VE + 3394\/2 + 4754, . . . 

(V2VE - V2) + (i (VE - 1) + 2) \/VE + 5 + V2VE + \/2, 4V2, 

2v / \ / 5 + 5 + iV2 (VE - l) , 2\/\/5 + 5 + i\/2 (\/5 - l) , 4\/2, 
i (V2VE - V2) + (i (VE - l) + 2) \/\/5 + 5 + \/2x/5 + \/2, 
[i (VE - l) + 4) \/\/5 + 5 + \/2\/5 + \/2, 
i (3V2\/5-\/2) + \/\/5 + 5 (VE + i (VE - l) +3) + 3V2VE + 7\/2, 
\/V5 + 5 (i (3\/5 + 5) + 6\/5 + 2) + i (8\/2\/5 + 8\/2) + 9\/2\/5 + 13\/2, 
\/\/5 + 5 (i (5\/5 + 5) + 12\/5 + 20) + 12 i\/2\/5 + 20 \/2\/5 + 46\/2, 
\/\/5 + 5(i (37\/5 + 95) +47\/5 + 65) + i (67\/2 VE + 187\/2~) + 81\/2\/5 + 137\/2, 
\/V5 + 5 (i (l68\/5 + 434) + 162 VE + 322) + i (332\/2\/5 + 796\/2) + 285\/2\/5 + 655\/2, 
\/\/5 + 5 (i (753\/5 + 1957) + 883\/5 + 1865) + i (l499\/2\/5 + 3579\/2) + 1644\/2\/5 + 3634\/2, .... 
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However, for other p on the unit circle, it is impossible to scale out denominators, even with a geometric progression. 
E.g., 



39i 27 r „ . , „ . „ r 39i 27 24i 47 
... i _ + _ >5> 3, + 4,3, + 4,5 > — + -, — + T> 

387 i 386 36783 i 3494 
25 + 25 ' 625 + 625 ' 
2088 i 48643 32642451 i 12616857 2038374144 i 1297298183 
25 + 625 ' 78125 78125 ' 1953125 1953125 

9109927539 i 22447118648 



1953125 9765625 



185 i 275 1or . 1or . ,„ ,„ 185 i 275 120 i 431 

. . . , 1 , 13, 5 « + 12, 5 1 + 12, 13, 1 , 1 , 

' 13 13 13 13 13 13 

5285 i 11722 4966665 i 4473518 

169 + 169 ' 28561 + 28561 ' 

468840 i 14119307 138123747005 i 79921772175 97867476702880 i 46503584483049 

2197 + 28561 ' 62748517 + 62748517 ' 10604499373 + 10604499373 ' 

438606726214525 i 3668262036619888 
+ 



10604499373 137858491849 ' 

wherein the powers of 1/5 and 1/13 grow quadratically. 

Interestingly, (lppl) permits definition of sequences containing 0. E.g., p := —1 gives a period 10 sequence with 
values in { — 1, 0, 1}. Alternatively, if p 3 + p 2 + 1 = p, the sequence is period 5. The appropriate generalization of 
Conjecture 4a is 



a k 



-i a k cifc+i cifc+2 (p 4 + 2p A + l) = (o|_! a 2 k+2 + a% a k+2 + a fc _i a A k+1 + a 2 k a 2 k+1 ) p 2 



There appears to be a generalization of this relation for initializations in violation of the palindrome property. 
Perhaps the most important of these is Sloane's A051138 

A-!,A , ... = -1, 0, 1, 1, -1, -5, -4, 29, 129, -65, . . . 

where 



A n — A— n — 



A n -\A n -3 + A 2 l _ 2 _ A n -i a n+ i — A n -2 a-n+2 _ A n _i a k -i a k -„ — A n -2 a- k a k -„-i 



We can think of A n as sinh and a n as cosh, but actually they're both theta functions. Also, a n is centered at 
n = 3/2 instead of 0. 

Solving this last equation for a k generalizes the Somos4 defining recurrence: 

A n+ l Cl fe _i a fe _ n _ 2 — A n+ 2 O-k-2 dfc-n-1 

a k = . 

Another such generalization is the "fc-tuple speedup" : 

A\ a n a n +4 k = A\ k a n+k a n+ z k — A k Az k a n+2 k ■ 
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Generalizing both of these is the three-variable relation 



A2 j Ak-\-j a n —j On— fc— 2j Aj -A-k+2 j 0jn—2j a n —k—j 



These expressions mixing A and a are somewhat striking because up until now, all the monomials in a given 
relation have had the same subscript sum, modulo A n — —A- n and a n — a 3 _ n . In particular, these nonconforming 
identities can not come directly from D s type determinant identities, except via the artifice of multiplying the 
deficient monomials by Ak and the overweight monomials by — A-k- 

A can be eliminated from the speedup identity via the relations 

A:$k _ a p dp+k+l a p+'.i fc+l fc ~ O-p+l Op+fc Op+3 fc dp+4 fc+l 

A k <lp+k a p +2 fc+l a P+3 k ttp+k+1 a p +2 k Op+3 fc+l 

. 2 2 2 

A 2 Am O-m + 2 fc + l a m+4 k ~ dm + 1 0- m +2 k a m+4 fc + l 

A\ Clm + fc a 2 m+ 2 fc + l a m+3 k ~ Om + fc + l ^^ + 2 fc a ™ + 3 fc + l ' 



for arbitrary m and p. Also, 

A k = at cifc+3 — afe+i cifc+2. 

If we eliminate A between this and the fc-tuple speedup identity, we get a polynomial in a with subscript sums 
which can be brought into agreement via selective application of On — 03 — n- 

Also, A n = S2n, where s n — Sloane's A006769: 

s- u s , ... = -1,0, 1, 1, -1, 1, 2, -1, -3, -5, 7, -4, -23, 29, ... , 

and 

Sn-lSn-3 + S„-2 
Sn — S — n — , 
Sn — 4 

the same recurrence as A„. Perhaps surprisingly, 

S2n+1 = (— l) n On + 2- 

That A1/2, A3/2, ■ ■ ■ can be integers suggests that 01/2,03/2,... could be, too. Substituting half- integers into the 
i9 expression below yields nonintegers, but it is likely that there are alternative analytic expressions for a n which 
disagree for nonintegers. 

Curiously, ^2 n does not obey the Somos4 recurrence. 

Note that a is even easier than A to eliminate from the mixed recurrences, since they hold for a = A\ I.e., 

Aj Ak A n A n _fc_3 j = A2 j Ak+j A n —j A n —k—2j Aj A^+2j A n —2j A n —k—j. 
With the relation A_ n = — A n along with linear changes of variable, this can be rewritten 

A—j Aj — k Aj — n A n + k+j ~t~ A2 j A — k A — n A n +k ~ A — j Ak+j A — n — k+j An+j, 

so that each term's subscript sum is 2j. We might thus expect an equivalent 3 by 3 determinant a la Conjecture 
4. The most general case gives a six variable relation with 24 terms of degree 6. The only apparent way to reduce 
to degree four is to specialize two of the variables to create terms of absolute value 1, i.e. A± 1 ,A± 2 , or A± 3 . 
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But this will introduce small integer offsets among the remaining subscripts, a feature notably absent from our 
trivariate relation. So for A n , at least, determinants may not tell the whole story. 

Likewise for #1: the trivariate relation empirically holds if we replace A n by #i(n, q), for arbitrary complex n and 
any fixed q within the unit circle. That it fails for the other $s suggests the existence of a four or more variable 
generalization. Indeed, by analogy with (4vars), 

Afc — i dk + i Aj — n dri+j ~ -Aj — i Clj + i Ak~ n a n -^-k -\- Afc—j Clk-\-j Ai — n tt n +j, 

also holding with A in place of a. So maybe (4vars) type determinants do tell the whole story. 
If s n = —s- n , si = 1, (as with s n := A n ,) then 



y, o, 

n, 0, 



-y s y+r 
-1 Sn+] 



s y — l 
-1 




Sn — y Sn-\-y S n Sy — \ Sy-\-l ~\~ S n —\ S n -\-l Sy 



(EDS) 



is equivalent to s being an elliptic divisibility sequence. Integer divisibility sequences merely require d\n =>■ Sd\s n , 
but the divisibility sequences discussed in this report appear to satisfy the stronger relation (s x ,s y ) — |s( Xiy )|, 
even when they disobey the addition formula. This may be what is meant by "strong divisibility sequence" . 

The EDS upside is this nice addition formula. 

Fomin and Zelevinsky have shown that Somos4, . . ., Somos7 are Laurent polynomials (rational functions with 
monomial denominators) in their initial values. 

Somos4oid polynomials: We can get true polynomials from the "odd" (s_ n = — s n ) sequences with the 
initialization — 1, 0, 1, 1, — 1, x, where x is unconstrained by the Somos4 recurrence, which gives 0/0. At greater 
length, 



s„ = —1,0, 1, 1, — 1, x,x + l,x 2 — x — 1, —x 3 — x — 1, — 3 a; 2 — 2 x,x° — x 4 + 3x 2 + 3 x + 1, 

-x 6 -2x 4 -5x a + 3x + l,~x 7 + 2x 6 -3x 5 -9x 4 -5x 3 -3x 2 -3x~l,.... 

A n is the case x = —5 and A n / 2 is the case x = 1. I. e., s„(— 5) = S2n(l). Empirically, this is a strong (redundant?) 
elliptic polynomial(l) division sequence for all x. If indeed the divisibility property holds for both integers and 
polynomials, then the values assumed by the polynomials Sk(x)/s^,n)(x) and s n (x)/s(k,n){x) are relatively prime 
for every integer x. 

It seems that the Chebychev polynomials U n -i(y) := sin(n arccos y)/ sin(arccos y) behave similarly. E.g., for 
integers k and n, sin(fcn arccos y)/ sin(n arccos y) is a polynomial in y, but of degree only (k — l)n. 

Here are the polynomial factorizations of s n (x) through n = 18. 
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n s n 

-1 -1 



1 1 

2 1 

3 -1 

4 x 

5 x + 1 

6 a; 2 - x - 1 

7 - (a; 3 + x+l) 

8 -a; (3a; + 2) 

9 x 5 -x 4 + 3x 2 + 3x + l 

10 -(x + 1) (x 5 -x 4 + 3x :i + 2x 2 -2x-l) 

11 - (a; 7 - 2x 6 + 3x 5 + 9x 4 + 5a; 3 + 3a; 2 + 3a: + l) 

12 -x (x 2 -x-l) (a; 6 + 2a; 4 + 5a; 3 + 9a; 2 + 9a: + 3) 

13 x 10 + x 9 + 13 x 7 + 12 x 6 - 6x 5 + 16a; 3 + 15 a; 2 + 6x + 1 

14 -(a; 3 +a.- + l) (a; 9 - 3 a; 8 + x 7 + 6 x 6 - 13 x 5 - 30 x 4 - 15 a; 3 + 4 a; 2 + 5 x + l) 

15 - (a,- + 1) (a; 13 - 2 a.- 12 + 5 a,- 11 - 2 a; 10 + 5x 9 + Sx 8 

-19a; 7 + 12a,- 6 + 63 x 5 + 50 x 4 + 20 a; 3 + 10 a; 2 + 5 x + l) 

16 x (3a; + 2) (2 a; 12 - 3 x 11 + 6 x 10 + 14 x 9 - 2 x 8 

+3x 7 + 23a; 6 + 18a; 5 - 6a; 4 - 27a; 3 - 27a; 2 - 12a; - 2) 

17 x 18 -3x 17 + 4 a; 16 + 6 a; 15 - 9 a; 14 + 5 a; 13 + 56 a; 12 + 69 a; 11 + 105 a; 10 + 311a,- 9 

+429 x 8 + 211 a; 7 - 2 a; 6 + 45 a; 5 + 135 a; 4 + 110 a; 3 + 45 a; 2 + 10 x + 1 



18 (a; 2 -a;-l) (a; 5 - x 4 + 3 a; 2 + 3 x + l) (x 1:i + x 12 + 7 x 11 + 19 x 10 + 25 x 9 

+7Sx 8 + 133a; 7 + 108a; 6 + 79 a; 5 + 65 a; 4 + 24 a; 3 - 6 a; 2 - 6 x - l) 

The degrees of the polynomials, starting with n — 1, go 

0, 0, 0, 1, 1, 2, 3, 2, 5, 6, 7, 9, 10, 12, 14, 14, 18, 20, 22, 25, 27, 30, 33, 34, 39, 42, 45, 49, 52, 56, 60, 62, 68, 72, 76, 81, ... , 

which is eight interlaced quadratic progressions: 

degs 89+ r = (4g + r)q + [-2, 0, 0, 0, 1, 1, 2, 3] r , < r < 7, 

which can be written 

degs„ = — (sin— )(sin— ) - (cos— )(- +cos— ) + — (2n - 5(-l) -5). 

It appears that n prime s n irreducible. The polynomials appear to be monic except for s 8n , whose leading 
coefficients appear to be (— ) n 3n. 

It appears that all the polynomials s n (x) have a root close to x — uj ^ —0.669499628215, + 21 uj + 

17 = 0, with proximity rapidly increasing with n. 

Besides the EDS condition, we retain the #i three-variable identity 

$2 j Sk S n Sn + fc = Sj Sfc — j Sn-j S n +fc+j "t - Sj Sfc+j Sn+j ^n+fc — j 
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This can be subscript-balanced as 

S2 j Sk+j S — n — k—j S n = S—j S — k Sj — n Sn+k+2 j S— j Sk-\-2 j S — n — k Sn+ j , 

but its asymmetry and failure to subsume the EDS condition suggest that we're missing a nice, /our-variable 
relation. Sure enough, by analogy with (4vars), 

Sk — i &k + i Sj — n Sn+j = Sj — i + i Sk — n ^n+fc &k—j Sk+j Si — n Sn+i 

withstands empirical testing. 

This identity specializes to a particularly attractive doubling formula: 

3 3 
S2n-1 = S n _ ± S n+1 — S„-2 S n 



S2r 



s n + 2 — S n -2 S^_|_i) S r , 



Given the four consecutive values s n ~2, • • • , s n +i, extend them to s n+ 3 stepping the recurrence twice. Then use 
the doubling formula to get the four values S2n-i, • • • , S2 n +2- Etc. 



We also have 



2 2 
Sb-a Sb+a — S a Sb-1 Sb+1 — Sa-1 «a+l «6 • 



And we have the ntuple speedup relation 

2 2 2 

Sfe S„ S n+ 4 k — S 2 k S n +k Sn+3 k ~ Sk S3 fc S n+2 k- 

This provides an alternative doubling process: Given four values s k , s 2 k, S3fc, s ik , start n at k and generate 
ssk, s&k, s-jk, s S k- Discard the odd multiples, and we have doubled k and are free to iterate. 

For x = 0, s n has period 8: 

sWO) = -1,0,1,1,-1,0,1,-1,-1,0,1,1,... 

= [0,l,(-l) 9 ,-l] r , 0<r<3, 

= &i~" 2/4 ??i(n7r/4,Q), « = 4q + r, b = .2653512762412i + .4652895036579, 

Q = .7359196601 139i + .3006597280279. 

Of course, a much simpler expression is 

s n (0) = sin (J^j - sin cos (J^j . 

For a; = — 1 the period is 5: 

amri-r(-l) = -1,0,1,1,-1,-1,0,1,1,... 

= [0,l,l,-l,-l] ri 0<r<4, 

= -^=i?i(2n7r/5, Q) n = 5q + r, b = 0.6155370356317, Q = .4856907848670i. 

vQ 
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For x = —2/3, we get eight interlaced progressions: 

s 8g+r (-2/3) = -1, 0, 1, 1, -1, -2/3, 1/3, 1/3 2 , -1/3 3 , 0, 1/3 5 , -1/3 6 , -1/3 7 , 2/3 9 , . 

[0, 3 3q+1 , (-)«3 2 « +1 , -3* +1 , {-) q+1 2, 3 1 -", (-)«3- 2 *- 1 , -3 1 " 39 ],. 



< r < 7, 



3 (2q+l)2 

= &M n2 i9i(ny,Q), n = 8q + r, b = -2.010659335767i, u = 0.8509811643954j, 
?/ = tt/2 - .7416161288587i, Q = .0026507066057. 

Note that t/ is not 7r/8 nor even real, so where do the periodic 0s come from? And \u\ is not 3 -1 ^ 16 , in fact, even 
its square root is too small. So where does the 3 ™ ' "growth" rate come from? The answer, as usual, is clear 
after Jacobi's imaginary transformation: 

2 g 9i 7T n 2 /16 ^ ^ttti^ 

s„(-2/3) = ^ 2 3„2 /16 8 . 

with g satisfying 

4i-2\ 1+j 
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^3 ' 
e.g., 

g « .5913080374704560258502159338438 + .4423170132359810537349781037012 i. 
We thus answer both questions, and reduce four mysterious parameters to one. Or rather, two, since 

q « .7241830710727415040344246937315 i + .5068861260317593704061905537186 
also satisfies the infinite product constraint, but produces the mysterious sequence 

VJ!,u,VdM, 3 , 3 ' 9 ' 9 ' 27' '243' 729' 6561' "" 



For x = the golden ratio, we get six interlaced progressions: 

S6q+r(<t>) = -1, 0, 1, 1, -1,0,0 2 ,O, ~<f> 4 , -(/>'', (j> 6 , -0 8 , -<^ 10 ,0, . . . 

= (-)V 6 (' , t 1 )[O,0- 29 ,0- 9 ,-l,0 < ' +1 ,0 29+2 ]„ 0<r<5. 

Likewise for the conjugate: 

s 6 ,+ r (-l/^) = -1,0, 1, 1, -1, -l/<t>, <t>~ 2 , 0, -<fi~ 4 , </>~ 5 , <t>~ 6 , -</>~ 8 , 0, . . . 

[0,<P«, {-</>)«, -l,(-4>)-<>-\ </>-*<>-% n<r< . 

(-)9^ 6 ( 2 ) 

= bu n2 -di(ny, Q), n = 6q + r, b = -1.403592671340i, u = 0.8476983265649z, 
y = tv/2 - .7507768082213i, Q = .0110573396552. 

This $ expression is close to s n (— 2/3) because —1/0 = —.618 is close to —2/3. 



13 



It is probable that s„(q) comes out in k such interlaced progressions when Sfc(a) = 0. E.g., when a 3 + a + 1 = 0, 
we appear to get seven interlaced progressions scaled by a ('=+ 1 H fe + :L /7) _ 

An algebraic x for which s„(x) has an elementary closed form is x = — w 3 — 2w 2 , where w 4 — w — 1, i.e. 
x 4 + 3 a; 3 - 5x 2 + 21 x + 17 = 0, w = -(3a; 3 - 26a; 2 - 89a; + 158)/283: 

s n = -1,0, 1, 1, -l,x,x + l,x 2 - x- 1,-x 3 -x-1, -x (3 a- + 2) , 17a- 3 - 38 a- 2 + 70 a; + 69, 
79 a; 3 - 126 a; 2 + 288 x + 273, -526 a; 3 + 925 a; 2 - 1838 x - 1803, . . . 
= -1, 0, 1, 1, -1, -w 3 - 2 w 2 , -to 3 - 2 w 2 + 1, 2 w 3 + 7 w 2 + 8 w + 3, 22 w 3 + 18w 2 + 25 w + 17, 
-w 3 - 11 w 2 - 24 w - 12, -465 u; 3 - 602 w 2 - 729 w - 389, 

-2073 w 3 - 2470 w 2 - 2983 w - 1653, 13809 u; 3 + 16717 w 2 + 20550 w + 11365, . . . 



(n-l) (n+1) 

sin I n arccos — - — I w 2 



, . , (n-l) (n + 1) 

f/n-l ^ W 2 



where U„ is the Chebychev polynomial, second kind. The degeneration of the •& corresponds to the vanishing of 
q. Note that one of the roots x ~ —0.669499628215, which is numerically close to —2/3, which at least explains 
the unusually small value of q in the otherwise puzzling $ expressions for s n (— 2/3) and s n (— 1/0). 

This #-free expression for s n (— w 3 — 2w 2 ) affords elementary expansions of b,u,y, and q about e = in 



s n (e -w 3 ~ 2w 2 ) = — u n ^ (ny, y/q) + 0{e 3 ), 



namely 



1 37_e _ __ 

1 + 757 w 3 + 450 w 2 + 87 w - 661 + 3241613 w 3 + 34180558 w 2 + 150075980 w + 30712993 + 
37 e 858967 e 2 

+ ■■■ 



31 x 3 + 86 x 2 + 192 x + 875 807948 a- 3 - 16234660 a; 2 - 29959955 x + 48768711 



Aw 2 



' ' ~ 16w 3 + 32to 2 + 2 ~ 1728448 w 3 + 518695360 w 2 + 706843904 w + 458996344 + ' 
= Vw I 1 + 



16 a- - 2 7072000 a; 3 + 1 10455488 a; 2 - 39784960 x + 144470392 

] 13j 

_ 308 ^ 3 + 200 u, 2 + 400^ + 258 ~^ + " 

y - arccos I 

13 e 290035289 e 2 



= arccos 



16 a- 3 - 28 x + 30 1143740480 a; 3 - 236851392 x 2 + 7982806560 a; + 9574893592 



2w 3 / 2 
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yy = 



q = 



qq 
qqq 



290035289 e 2 



32238208032 to 3 + 35449717760 w 2 + 30488323136 w + 11959840616 ' 

/ 4339 £ 

6 V 1733 w 3 - 55408 w 2 + 2992 w + 36867 +qq + ' 
234- 77w + 91w 2 -170w 3 

/ 4339 e 

6 V + 1696 x :i + 4926 x 2 - 28957 x - 26043 ~ qqq + ' 



13 x 3 + 31 x 2 - 72 x + 344 



15857700151 to 3 - 1431606991 w 2 - 2016123508 w + 7704208617 ' 



915645540 a; 3 + 3113407751 x 2 - 257448438 x - 3676219901 ' 



What is it with 283 = 566/2 = V80089 = \/22665187? Answer: -283 = discriminant (w 4 - w - 1) = 
^/discriminant (x 4 + 3 x 3 - 5 x 2 + 21 x + 17) . 

Using the negative root w w —0.72449195900052, these expansions through e 4 go 

6 w -0.9125730603509 i (1 - 0.04769878803144 e- 0.01401078082448 e 2 

-0.0050002222638 e 3 - 7.44689129426338 • 10~ 4 e 4 H ), 

u w 0.85117093406702 i(l- 0.07866586437604 e- 0.01169998061242 e 2 

-0.00163136901344 e 3 + 1.1997262491263 • 10 -4 e 4 + • • •), 
y w arccos(0.81081103497608 i (1 + 0.29582732047992 e + 0.07670823278512 e 2 

+0.02128110984092 e 3 + 0.00440366306328 e 4 + • ■ ■)), 
g w 0.00248633800734e(l - 0.87527986918762 e + 0.20054769983458 e 2 + 0.0025905870597 e 3 + ■ ■ ■), 



suggesting a fairly commodious radius of convergence. 

The positive root w = 1.22074408460576, s n (e - w 3 - 2w 2 ) = -1, 0, 1, 1, -1, e - 4.7996, . . . , has expansions 

6 = .97451669219348 (1 + . 02478371881502 e+. 00317253690858 e 2 

+4.97976424690216 • 10~ 4 e 3 + 8.74003233413136 • 10" 5 e 4 + ■■■), 
u = 1.10487288165008 (1 - .01269137385174 e - .00121848904668 e 

-1.63711992744874 • 10~ 4 e 3 - 2.5512354068948 • 10" 5 e 4 + ■ ■ ■), 
y = arccos (.37070894172584 (1 - .00810147538142 e - . 0018053537704 e 2 

-3.84723167907628 • 10 -4 e 3 - 8.0356533846957 • 10" 5 e 4 + ••■)) , 
q = -.02972037154846 e(l + . 1115487179065 e+. 01495508864096 e 

+ .00221611060456 e 3 + •••)• 

Another interesting Somos4 (apparently (strong) polynomial (non-E)DS) is 

so, si, ■ ■ ■ = 0, 1, i, 1, x, i (x — i) , — i (a; 2 + ix + l) ,—i (a; 3 — x + i) , — i x (3 x — 2i) , i 
(x 5 + i a; 4 + 3 i x 2 + 3 x — tj , i (x — i) (a; 5 + ix 4 — 3x 3 + 2i a; 2 — 2a; + i) , 
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—i [x 7 + 2i x 6 — 3x J + 9ix 4 + 5 x 3 ~ Six 2 — 3x + i^j , 
—x (x 2 + i x + l) (a; 6 — 2 x 4 + 5 i x 3 + 9 x 2 — 9 i x — 3) , 

which gives us Gaussian integers, among other things. As with the previous Somos4 polynomial sequence, there 
is likely a value of x for which the #i degenerates to a Chebychev, and consequently another set of elementary 
expansions of the i9 parameters about this x. But foo, these polynomials are essentially identical to those generated 
by the —1, 0, 1, 1, —1, x,x + 1, . . . sequence. 

Corollary 4: the determinant 

D a ( S > *' U > V )=0, 
\w, x, y, z J 

where s,t,u,v,w,x,y, and z are arbitrary integers. Proof: Dodgson's rule, provided the central 2 by 2 doesn't 
vanish. 

Expression as Email from Noam Elkies to sci.math suggests the relation 

oo 
k— — oc 

Using n € {2, 3, 4, 5} to numerically approximate b, u, q, and z, 

{b = 1.01943271913292, u = 0.63853138366726, z = 0.05462469648874, q = 0.02157360406362} . 
These constants do not appear to be in Plouffe's collection. Plugging in —3, —2, . . . , 10, 11 for n gives 
{6.99999999999998, 3.0, 2.0, 1.0, 1.0, 1.0, 1.0, 2.0, 3.0, 6.99999999999998, 23.0, 58.9999999999998, 

313.999999999998, 1528.99999999998, 8208.99999999994}, 

in good agreement. Due to the pleasantly small value of q, we even get 1123424582770.98 for ars = 1123424582771. 
In fact, the only terms affecting this double precision result were —9 < n < —2, making the series a competitive 
alternative numerical method. 

With Jacobi's imaginary transformation, we get an even nicer, entirely real expression: 

a n = bu (n ~ 3/2) ^ 4 (( n -^) V'l) 

with 

y w 1.9511889024071, q w .07632928490026, 

# 4 (ZJL q ) 1/8 I ^ 4 (M. „) 
b = — V 2 - ; o/c w .92252487906093, u = J ) 2 ' { w 1.10763024250632. 

Mu) 1 V^(^«) 

A n not surprisingly comes out as a #1, also with a fairly small q: 

oo 

A n = -2bu n2 ^](-l) fc q (2k+1)2 sin((2fc + l) ny) 

k=0 

= -bu" &i (ny,q 4 ) , 

where 

y = 1.9511889024071, § = 0.52562110924304, b = .92252487906093, u = -1.10763024250632. 
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Note that this q is raised to the fourth power in the $1, so that, modulo an alternating sign, these parameters 
are identical to those in the $4 formula for a n . I found these parameters enormously tough to compute (prior 
to Jacobi-transforming the a n expression), which may explain their absence from Elkies's email. Then again 
I flunked numerical analysis. Of course, now that we have it, the Chebychev expansion is also valid at w = 
1.22074408460576, x = w z - 2w 2 = -4.79960475359606, which is close enough to -5 to provide an excellent first 
approximation. To a different solution, however! (Negated b and u, 7r-complement of y.) 



-4.99999999999998d0, -4 = -3.99999999999998d0, 29 
-65.0000000000002d0, -3689 = -3688.99999999996d0 



Testing the non-Chebychev expression: 

= O.OdO, 1 = l.OdO, 1 = l.OdO, -1 = -l.OdO, -5 = 
28.9999999999998d0, 129 = 128.999999999998d0, -65 : 

Recalling that A n — s n (— 5) = S2 n (l), where s„{x) is the EDS polynomial sequence satisfying Somos4, we sought 
a $1 expression for s„(l) to see if S2 n (l) gives the same expression as the #1 for A n . In fact, we found (with 
much difficulty) ten $1 expressions that agree with s„(l) for integer n. For noninteger or nonreal n, symmetry 
suggests that there are as many as sixteen different functions. (And not one of them coincides, for n <— 2n, 
with our A n expression.) The sixteen seem to divide into two classes of eight. Within each class, their values at 
n = 1/2, 3/2, . . . agree modulo conjugation and multiplication by some integer power of i. 

However, there appear to be more than sixteen (y, q) pairs producing s„(l). I.e., there are multiple ways to 
express s„(l) as bu n •d\{ny,q) that agree even for complex n! In particular, 



b = 



tfr (y, q y #i(3y,q) 



$i(2y,q) 3 ' 



and 



g $i(2y,q) 3 Vi(ny,q) ( gi (y, gf gi (3 y, q) 

&i(y, q f #i{3y,q) V $i(2y,qf 

seems to be exactly the same function of complex n for 

y w 0.49235539271999 - 0.74875275029651 i, q^ 0.69018582634555 i + 0.54229640598463 

as for 

y w 1.11453161008963 i + 0.62943384983216, q w 0.63418111840451 i - 0.43035475675355. 
(These are not mutual Jacobi transformations.) 

Over a period of several hours, an automated grid search cum Newton's method turned up the following approx- 
imate (y,q) pairs for s n (l): 



V 



q 



1.2482046102601 i + 0.04657952537373 
0.78322226431624 i + 0.30719109513916 
0.49235539271999 - 0.74875275029651 i 
1.0429573809123 i + 1.37265463822724 
1.00790006282379 i + 0.75028009770783 
1.11453161008963 i + 0.62943384983216 
0.2917595098002 i + 1.21452080661459 



0.2041895179564 i + 0.06533908137423 
0.45270853094805 i + 0.20573877584467 
0.69018582634555 i + 0.54229640598463 
0.73581195373709 i + 0.23968436737044 
0.73581195373709 i - 0.23968436737042 
0.63418111840451 i - 0.43035475675355 
0.2041895179564 i + 0.06533908137423 



0.4615909499519 i + 1.463034339711 0.86308985589011 i + 0.319843560314 
2.01129861135631 i + 1.12816860769464 0.45270853094804 i + 0.20573877584469 
0.23882852142275 i + 1.76522302148841 0.73581195373708 i + 0.23968436737043 . 
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The last three are quite unlike the first seven for noninteger n. For integer n, s n (l) is real and the conjugates 
of all these work as well. Note the equality of the first and last q in the first group. This would seem to be an 
instance of translation by a quasiperiod, via 

„ / • , n (-!)" e 2tny i) 1 (y,q) 

01 (y + tnlogq,q) = ^ -g 

q n 

for some integer n. But there are three problems with this. First, following a tradition that still puzzles me, 
we have made no provision for a geometric (r n ) factor in our $ formulae, even though the recurrence relation is 
unaffected by such a factor. (Ah, but the EDS relation is affected.) Second, if we go ahead and solve for n, 

n = 0.75842109957414 i. 

Not an integer, but pure imaginary, for some reason. Third, if, for some integer n, one of these turned up in our 
search, why wouldn't we find dozens more engendered by other values of n? 

Also compare the fourth and fifth q of the first group with the last q of the second group. This offers hope for 
some simple relation between the corresponding y. 

The grid search also turned up three spurious pairs, 

v q 

0.16973145507896 i + 0.59439464562658 0.7091459745365 i + 0.51937050102005 
0.15670471616132 i + 0.57764641771107 0.61450698230835 i - 0.63841374520714 
0.13574561978509 i + 3.27166894673346 0.39191084430236 i + 0.91218127829481, 

which generate the sequence 

_ 144r n - 3 r n -i +432r^_ 2 

Tji — V— ji — 

= . . . , -1, 0, 1, 2 2 3, -2* 3*, 2 7 3 6 , 2 12 3 10 , -2 16 3 15 , -2 23 3 20 , -2 29 3 26 5, 2 37 3 33 7, 2 46 3 41 , -2 55 3 50 13, 

2 6S 3 61 ^ 2 77 3 70 2 90 38! ^ _ 2 103 g 93 ^ _£VI 3IO6 g 53^ 2 133 3I2O n ^ _ 2 148 gl35 ? m ^ 



whose significance thus far eludes me, although it appears to be a (weak) EDS. 

Somos5 differs from Somos4 in two respects: b n — b^-n instead of a n = a3_ n , and a different order-reducing 
relation from Conjecture 4a. 



Substituting into Conjecture 4.5 s„ = b n ,u = x,y = 0, z 

, a; + 1/2, 3/2, 1/2 
1 x- 1/2, -1/2, -3/2 



-1, v = 1, w = 0, 



bi b*2 x 


b x b x +i 


b x -i bx+2 


2-x bx + 1 


bi b 2 


b b 3 


bi-x b x 


bo bi 


b-i b 2 


b2x 


b x bx+i 


bx-i bx+2 


x + 1 bx+2 


1 


1 


b x bx+3 


1 


2 



= 0. 



This gives us 62^ in terms of five values near b x - Alternatively, put a; — 1 for a;, then u = x,y = 2,z = 0,v = 0,w = l 
to get 



D h 



x + 1/2, 1/2, 3/2 
a; -3/2, 3/2, -1/2 



&2&2a;-l b x -l b x + 2 b x b x + l 

b 2 -x bx-i b-i b 2 bo bi 
b-A-x bx bo bi bi 62 
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&2x-l 

b x -i b x+ 2 
b x b x +i 



-1 bx+2 

2 
1 



b x b x +i 
1 
1 



= 0. 



This gives us b2 X -i in terms of four consecutive b values. However, we can reduce these to threel 



D b 



fc-1/2, fe + 1/2, fc + 3/2\ 
fc-1/2, fc + 1/2, fc-3/2 J 



&2fc-l 2&2fe 62 fc-2 
&2 fe 62 fc+1 62 fc-1 
&2fe+l 6 2 fc+2 62 fc 



= 0, 



, &2 fc-2 &2fe+l ~~ 3&2fe 62 fc-1 62 fc+1 + 2 fojfc 

»2fc+2 — To 7 7 > 

°2 fc-1 — °2 fc 02 fc-2 

a fourth order recurrence. Alternatively, 

/ fc-3/2, fc-1/2, fc + l/2\ 
6 ^fc-1/2, fe + 1/2, k + 3/2j 

— 2 ((3 62 fc-1 — 4 62 fc 62 fc-2) 6 2 fc+2 + 3 6 2 fc-2 &2 fc+1 — 76 2 fc 62 fc-1 62 fc+1 + 5&2fc) 

or 

_ 3 62 fc-2 &2fc+i — 7&2fc 62 fc-i 62 fc+i + 5 6 2fc 

2fc+2 3fe| fc _ 1 — 4&2fc &2 fc-2 

a different fourth order recurrence. Subtracting, 

&2 fc-2 &2 fc + 1 + (2&2fc-l — 5&2fc 6 2 fc-2 &2k-l) 62 fc + 1 — b\ k b\ k _ 1 + 3 62 fc 62 fc-2, 

a third order recurrence for &2fc+i in terms of the three previous terms. (Assuming you know which sign to take 
on the square root). But what about &2fc? Simply replace fc by 3 — fc and b x by &4_x, and we have &2fc in terms of 
b2fc-i,&2fc-2, and &2fc-3: 

3 2 2 3 2 2 

2£>2fc-3 62 fc-1 — &2 fc-2 b 2 fc-1 — 5 &2 fc &2fc-3 &2 fc-2 &2fc-l + 3 62 fc &2fc-2 + &2 fc &2fc-3 = 

But this sensitivity mod 2 entails four residue classes when we order-reduce the duplication formulae.: 

22 33 222 4 

= b 2x+ 2 &4x-i + 6 b 2x b 2x +2 &4x-i — 21 6 2 x 6 2 x+i b2 X +2 &4x-i + 16 62a; b 2x+1 b 2x+2 &4x-i 
-46 

x+i bix-i + 9 &2 x b 2 x+2 — 36 b 2 x b 2x+1 b 2 x+2 + 57 6 2 x &2x+i b 2x+2 — 36 b 2 x b 2 x+1 b 2 x+2 
+8b 2x b 2x+1 

= — 3 b\ x b 2 x+2 + 4 b 2 x b 2 x+1 b 2x+2 — 2b 2x+1 + 64 x 

22 33 222 4 

= b 2x b 4x+1 + 6 b 2x b 2x+2 64 x+i — 21 b 2x b 2x+1 b 2x+2 b4 X +i + 16 b 2x b 2x+1 b 2x+2 641+1 

-4 62 x+1 64 x+1 +9 62a: 62 x+2 — 36 62i &2 x+1 62 x+2 + 57 62 x bix+l 621+2 — 36 62 a: 62 x + 1 62 x+2 
+8 62 x + 1 62 x+2 

= b\ x 64 x+2 — 6 &2 x 62 x+2 64 x+2 + 32 62 x b\ x+1 62 x+2 64 x+2 — 62 62 x b 2 x+1 b\ x+2 64 x+2 

+40 62 x 62 x + i 62 x+2 64 x+2 — 8 &! x+1 64 x+2 + 9 &2 x &2 x+2 — 48 62 x b\ x+1 b\ x+2 
+86 62 x 62 x + 1 62 x + 2 — 56 62 x 62 x + 1 62 x + 2 + 12 62 x + 1 62 x + 2- 



We can obviate the first or last of these with (respectively) the odd or even version of the third order recurrence. 
This takes care of doubling. 
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As with Somos4, we assume three or four values near b„ x and another tuple near 6(„ + i) x . Then 



D h 



(n + l)x + 1/2, 
nx - 1/2, 



1/2, 
-3/2, 



3/2 
-1/2 



bx + lb(2n+l)x &(n+l) x-l &(n+l) ic+2 &(n+l) a: fyn+l) m+1 

bnxb n x+3 2 1 

bnx+1 b nx +2 1 1 







gives &(2n+i)a;- Similar constructions provide the adjacent values, and in principle, we can use the third order 
relations to make everything work on triples. 

Note, however, that we could avoid the square roots and mod 4 intricacies by maintaining four values instead of 
three, with the help of the fourth order relations that we subtracted to get the third order one. 

(Brief flame: a nearly forgotten fact of hardware design is that a binary square root instruction via the "schoolboy 
algorithm" is actually simpler than the divide instruction. In the early 1960s, the Packard-Bell 250, as feeble a 
machine as you could imagine, whose active registers were magnetostrictive delay lines instead of flip-flop words, 
and whose divide instruction needed a software followup correction, nevertheless had a hardware square root (with 
remainder) that worked perfectly, in the same time as an uncorrected divide.) 

But which sign of the square root do we take? Not obvious! E.g., suppose we try to use the third order Somos4 
relation to compute a x from the three previous values: 



\/-^ a l-3 a l-l + l2a l-3 a l-2 a l-l -8^-3^-2^-1+^-2 . 2^-2^-1 



+ ■ 



3 

a x-2 



2<_ 3 Clx-3 2o^_ 3 

Then for 2 < x < 38, the sign s x coincides with 

sgn ((f) (x - 1) - rnd (<f> (x - 1))) , 

where <fr is the golden ratio, and rnd(x) := [x + 1/2J, the "round" function. But for x = 39, this fails! 

In practice this isn't really a problem, since we can simply choose whichever produces an integer value for a x , and 
we can usually check this modulo something small. But there may be another solution. 

To take a simpler example, from the third order relation for b„, we notice that 



v/ fc 2n+i&2n+3 -&ln+2 = °> h 8 > 57 > 455 > 22352, 47767, 69739671, . . . 
may be an elliptic divisibity sequence. Trying various sign patterns, we eventually find the recurrence 

h„ = _ 8fe "-4^-i+57fo„- 3 fo„- 2 = ^ Q ^ _^ ^ g ^ 5 ^ _ 45 ^ 22352 ^ 477 ^ 69739671) 33 85862936i .... 

h n -5 

But this disobeys the EDS formula. Searching further, 

g„ = 57 9n - 39n ~ 2 -8gn-4ffn-i = _ 1; 0j 1; _ 1; _ 8; 57; 455i __22352, -47767, 69739671, -3385862936, . . . 

g-n-5 

appears to satisfy the EDS addition formula. Evidently, g n — —h n except when n — (0mod4). 

It seems reasonable to conjecture that the desired sign pattern for these square root expressions is one that yields 
an EDS, or at least a simple recurrence. An example of the latter is 
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v / 3b 2n b 2n+2 -2b 2 2n+1 = 1, 1, 1, 1, 7, 1, 391, 2729, 175111, 
which cannot be an EDS. But there is an assignation of signs: 

/„ = 57 /"- 3 /"-2 ~8/n-4/u-i = _ lj _^ ^ ^ _^ _^ 391] _ 2729; _ 175111] 8888873, 565353361, . . . , 

Jn-5 

i.e., the same recurrence as <?„, above. 

Unfortunately, the order of the recurrence is likely to exceed the order of the relation that engendered the square 
root, defeating the presumable purpose of maintaining smaller intervals of consecutive values. 

However, these / and g sequences serve another purpose if we interlace them: 

D Bn-4 B n -\ + B n -3 B n -2 D 

Dn — r= — —O-n 

-Dn-5 

= -1, 0, 1, 1, 1, -1, -7, -8, -1, 57, 391, 455, -2729, -22352, -175111, -47767, 8888873, 



which is not quite an EDS. Yet we have 

k _ B2n+1 fcfc-1 bk-2n-2 — -E?2n+2 frfc-2 &fc-2n-l 
B2 n bk-2n-3 

generalizing the Somos5 recurrence. 

Note that B n obeys the Somos5 recurrence, yet somehow jumps from 1 to 7, via the well known identity —7 = 0/0. 
(Of course, in the limit of absurdity, any sequence which is alternately satisfies Somos5 and Somos7.) 

The fc-tuple speedup formula is 

Bh £?2 k b n b n +5 k = B2 k Bg k b n +k b n +4 k — Bk B4 k b n +2 k b n +3 k ■ 

A three-variable generalization of these last two relations is 

^ _ Blj B2k+j b n -j b n -2k-2j ~ Bj B2 k + 2 j &ra-2j b n -2k-j 

Bj B 2 k b n -2k-3j 

This is the same relation as with a and A, except that k must be even. The further generalization k <— fc/2 
involves an additional term that appears guessable, and seems to vanish for even j. 

Thus, if b n is the Somos5 cosh, then B n is the sinh, although a glitch is that b n is centered at n — 2 rather than 
n = 0. As with a and A, we can merely substitute B for b in the last identity: 

Bj B 2 k B n B n -2 fc-3 j — B2 j B2 k+j Bn-j B n -2 fe-2 j ~ Bj B2 k+2 j B n -2j B n -2 k-j- 

This all suggests that A and B will have cheaper addition algorithms than a and b (Somos4 and 5). But caution: 
even though B„ isn't quite an EDS, its values mod 19 lie in {0, 1, 7, 8, 11, 12, 18}. There may be other moduli with 
even (proportionately) sparser residue classes. 

As with A4, we can replace the B5 — —7 (= 0/0) term by x to get a sequence of polynomials: 

s n = — 1,0, 1, 1,1,— l,x, x — 1, — l,x 2 — x + 1, —x 3 + x 2 — 1, —x (x 2 — 2x + 2) , —a; 4 + x A — 2x + 1, 

(x - 1) (x A - x 3 + x 2 + l) , -x 6 + 3 a; 5 - 3 a: 4 + 3 a; 2 - 2x + 1,2 a; 5 - 5 a; 4 + 6 a: 3 - 2 a; 2 - x + 1, . . . . 
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These do not form an EDS, even with reassignation of signs, except when x = —2. But they retain "strong" 
divisibility, 

(s k (x),s n (x)) = |s(fc,„)(a;)|, 

for all x. 

What does (apparently) form a (weak) EDS is the (algebraic) sequence 

(_l)n + n 2 

t„ := s„ (-X - 1) 8 , 
which could be made (peculiar) polynomials with x = — 1 — y 8 : 

-, n 1 5 8 17 24 / 8 , ,\ 37 / 8 . \ 48 65 I 16 , 8 , \ 

■■■,-1,0, l,y ,y ,-y ,-y (y + 1) , -y (y +2J ,-y ,y (y + 3y + 3), 

80 / 24 . , 16 . r 8 . , \ 101 I 8 . -, \ / 16 . , 8 . r\ 

V (y +4y +5y +l),y [y + 1) (y +4y +5), 

120 / 32 . r 24 . n 16 i r 8 -, \ 

-j/ (y +5y +9y +5y -1) 



This EDS satisfies both the Somos5oid 

tn-3 tn-2 {—X - 1) 3/2 + tn-A tn-1 (-S - 1) 



tn-5 

and the Somos4oid 

_ tn-3 tn-1 (-» ~ 1) 5/4 ~ fn-2 (-» ~ 1) 
*n — 

tn — A 

Eliminating t n _ 5 and t n -4, 

(t a n -2 tn + t n -z in-l) ( — X — 1)'' /4 + t n -3 t n -2 t n -l t n \/ — X — 1 + t n _ 3 t% 
= *n-3 tn-2 tn-l t„ ( — X — 1) 3/2 + t„- 2 t„_i ( — X — 1) . 

Solving for t„ yields a radical whose sign seems to depend on x. 

In the special case x 3 + 5x 2 — Wx + 11 = 0, this EDS has the elementary formula 



,'41-35a:-x 2 \ T ~ TT j 1 s/llx 2 + 17a;- 244 

tn = — T t/„-l 



23 / I 2 V 23 / 

This is basis for the Chebychev expansion for the "sinh" analog of Somos5. 

As in the Somos4 analog, the s pr imo polynomials are irreducible, at least through sqi. 

The polynomial degrees go 

0, 0, 0, 0, 1, 1, 0, 2, 3, 3, 4, 5, 6, 5, 8, 9, 10, 11, 13, 14, 14, 17, 19, 20, 22, 24, 26, 26, ... , 
being fourteen interlaced quadratic progressions: 

deg s 14q+r = (7q + r)q + [-2, 0, 0, 0, 0, 1, 1, 0, 2, 3, 3, 4, 5, 6] r , < r < 13. 
The polynomials appear to be monic, except for S7n, whose leading coefficients appear to be ±n. 
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Happily, the three-variable relation seems to hold for general x: 

Sj S'2k S„S n _2fc-3j = S2j S'2k+j S„-j S n -2k-2j ~ Sj S2 k+2 j S„-2 j S„-2k-j- 



Better yet, putting u = —x, v = —y and s„ = — s_ n in Conjecture 4.5, we get the four (integer) variable, three 
term relation 

&k— ;7+l ^fc+j S— n— S n — i = S— j — Sj — i S — n +k + l ^n+fc T" -Sfc — i+1 S — n — j + 1 Sn—j • 

Caution: This fails for noninteger i,k,j,n even though the subscripts are integral. 

In his email to sci.math, Elkies makes the remarkable observation (modulo typos) that t n := (2/3)'" mod 2 ^ 4 b„ 
satisfies the reduced order (quasiSomos4) recurrence 

t„-2t n +2 = V6t n -!t n+1 - t 2 n (Elkies). 

It is probable that t n also satisfies a third order relation (of higher degree). 

As with the Somos4 polynomials, it appears that s n (a) falls into k interlaced elementary progressions when 
Sfe(a) = 0, but they are more complicated. Alternatively, they can be written as mk simpler progressions, for 
some multiple m. E.g., ss(e™/ 3 ) = and s n (e" r / 3 ) is merely periodic, but the period is forty-eight! 

Also like the Somos4 polynomials, we can get a #i expression via the change of variables: 

tan (arctan ty-x - 1 + ^f) 

tn(x) — — = S n (x), 

\/ — X — 1 

where t n (x) satisfies 

tn-4. tn = t n -3 tn-1 V~X — 1 — t^-2- 

It shouldn't be hard to find a Chebychev formula for some algebraic x, and hence elementary expansions for the 
# parameters, as we did for Somos4. 

Expression as Elkies' email gives 

oc 

&„ = (3/2) (nmod2)/ W"- 2)2 J2 1 k2 z k(n ' 2 \ 

k — — oo 

with q = 0.02208942811097933557356088 . . . , z = 0.1141942041600238048921321 

b = 0.9576898995913810138013844 . . . , and u = 0.7889128685374661530379575 .... These constants do not appear 
to be in Plouffe's collection. 

Similarly, 

s„ = B n tan (j-^- + arctan 
also satisfies the (Elkies) recurrence, giving 



Bn = b cot arctan u n ^(-lf q {2k+1) sin ((2 k + 1) ny) 

fc=0 

— — cot (j-^- + arctan Voj u n #i (ny,q A ^ , 
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where 

6 = -1.82905778669392, it = -1.07425451486466, y = 0.89396990235568, q = -0.52353014451686. 



Numerically testing this equation for — 1 < n < 22: 

-1 = -0.99999999999986d0, = O.OdO, 1 = 0.99999999999986d0, 1 = 0.99999999999986d0, 

1 = 0.99999999999984d0, -1 = -0.99999999999984d0, -7 = -6.99999999999914d0, -8 = -7.99999999999886d0, 

- 1 = -0.99999999999982d0, 57 = 56.9999999999918rf0, 391 = 390.999999999952d0, 455 = 454.999999999928d0, 

- 2729 = -2728.99999999958d0, -22352 = -22351. 9999999972d0, -175111 = -175110. 999999977d0, 

- 47767 = -47766.9999999918d0, 8888873 = 8888872.99999885d0, 69739671 = 6.9739670999992d + 7, 
565353361 = 5.65353360999916d + 8, -3385862936 = -3.3858629359995d + 9, 

- 195345149609 = -1.9534514960898d + 11, -1747973613295 = -1. 74797361329478c? + 12, 

- 4686154246801 = -4.6861542468004d + 12, 632038062613231 = 6.32038062613152d + 14. 

In email to sci.math, Randall Rathbun and Ralph Buchholz make the remarkable claim that the Heron triangles 
with two rational medians have side lengths 

[Bi+3 (Bf bf +3 bj +4 + Bf +1 bf +2 Bf +2 ) h+5, 
B i+ 2 b i+ 4 (B 2 bf +2 Bf +3 6 2 +5 + B 2 +1 B 2 +2 b 2 + 3 b 2 +4 ^ , 

Somos6: 

2 

^n— 3 Cn— 4 C n — 2 Cn— 5 C n — \ C n — Q Cn 

13 5 9 

11 3 5 

11 13 

= — 4c n _ 6 C n + 4c„_ 5 C n _l + 4c n _4 C„-2 + 4c 2 _ 3 , 



/n-3, 0, 1, 2\ 
c ^ 0, 1, 2, 3 J - 



(four times) the defining recurrence for Somos6. But 



A 



= 80 = A 



0, 2, 4, 6 
0, 1, 3, 4 

so 4 by 4 isn't enough. Building on the nonsingular matrix, 



A 



x, 0, 2, 

y, o, l, 



4, 6 
3, 4 



Cx — y Cy+x 

C-y Cy 

C'2-y C y + 2 
C-A-y Cy+4 

CQ — y Cy + ft 



_ I 

1 2 

2 ' 



-1 Cx+1 

3 
1 
1 

5 



Cx-3 C x +3 
9 

3 
5 

23 



-4 C^+4 

23 
15 
9 
75 



= 80 C x -y C y+X + (4 C x -4, C x +4 + 12 Cx-3 C x +3 ~ 52 Cz-i Cx + 1 — 44 c 2 ) C y _i C H+ 6 
+ (— 16&E-4 Ca; + 4 — 88Ca;_3 C x +3 + 328 C x -1 C x + 1 + 176 C 2 ) Cy + i C y+ 4 

+ (-28 c^-4 c^+4 - 44 Cx-3 Cx+3 + 284 c x -i c x +i + 188 c 2 ) c H+2 

+ (8^-4 Cx+4 + 2Ac x -3 Cx+3 — I44c x -1 Cx + 1 — 48 c 2 ) C y +5 Cy, 
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which empirically vanishes, providing a fairly messy addition formula. Using the defining recurrence to alge- 
braically eliminate c y +a, c x -4, etc., yields even messier relations, but with lower order (= width). The determinant 
approach can't do much better unless we can find a "narrower" non-singular 4 by 4. But this is unlikely, since re- 
placing the sixteen numerical coefficients by undetermined ones indicates that the c y -iCy+e term is indispensable, 
and the c x -4C x +4 term is not replaceable by Cx-2Cx+2. Slightly better may be 



80 C x - y +\ C y + X = ( — 4 C x -3 C x +4 — 8 C x -2 C x +3 + f 6 C x -i C x +2 + 28 C x Cx+l) Cy-2 Cy+H 

+ (32 c x -3 c x+ 4 + 24 c x -2 c x+ 3 — 88 c x -i c x+ 2 — f 44 c x c x +i) c y -i c y+5 
+ (-44 c x - 3 c x+ 4 — 48 c x -2 c x+ 3 + 176 c x -i c x+2 + 1 88 c x c x+ i) c y c y+i 

+ (8 C x -3 C x+ 4 + 96 C x -2 C x+ 3 — 1 52 C x -l C x+ 2 — 96 C x C x+ i) C y+ i C y+ 3 

from 

fx+ - - 1 

D f ~ 2 ' 2' 2' 2' 2|_n 

\ It 2' 2' 2' 2' 2/ 

Conjecture: there is no sequence of bivariate polynomials obeying Somos6. Evidence: Initializing with 

n t 3 6252 8 3 I 4 , -,\ 15 5 

s , si,... = 0, l,x y, -x y ,x y ,-x y [xy + l),-x y,..., 

gives polynomials through S29, but then S30 has a denominator of x. Replacing y <— yx will probably move the 
violation a few terms to the right. 

However, if we initialize with 

s ,si,... =0,1,1,1,1,-2,0;,..., 

then we appear to get polynomials 

s 7 ,s 8 , . . . = x - 1, 2x + 3, x 2 + 5, x 2 + x - 9, - (x a + 2x 2 + 4 x + 2) , x 2 + 13 x - 13, 
- (x 4 + 2x :i + 7x 2 + 7x + 32) ,x 4 - x 3 + 10 x 2 -41a;- 13,x 5 + 5o; 4 + 10 x 3 + 7 x 2 + 32 x + 70, . . . 

which is decidedly not a divisibility sequence. In fact, the only reducible s n through sqq is 

sie = -(0; + 1) (3 a; 4 -x 3 + 21 x 2 - 52 x + 117) ! 

And, unlike Somos4 and 5, s fe (a) = 7^ s mk (a) = 0. 
Somos7: Random clue: 

dfc-11 d k + dk-10 dk-i + d fc _7 d fc _4 = 61 d fe _6 rffc-5- 



D d 



1 

2 ' 



1 

2 ' 



3 5 

V ? 

2 ' 2 



dn — 4 dn-3 


d n -h d n — 2 


dn — Q dn — 1 


d n —7 dn 


1 


3 


5 


9 


1 


1 


3 


5 


1 


1 


1 


3 


-4 (d„_7 d n 


■ dn — 1 


dn — 5 d n - 


-2 — d n _4 



dn-3) , 



the defining recurrence. But 



Ai ( _ I ' I ' r ) wo, 

2 ' 2 ' 2 ' 
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so 



d x -id x +2 d x d x +i d x -2d d x +5 

9 15 17 137 

3 1 5 17 

1 1 1 15 

113 9 

and 

160 d x - y +\ d y + x + (4 d x -4 d x +s + 12 d x -2 d x +z — 28 d x -i d x +2 — 44 d x d x +i) d y -2 d y +7 
+ (12 d x -A d x +5 + 36 d x -2 d x +z — 164 d x -\ d x +2 — 52 d x d x +i) d y d y +$ 
+ (—28 d x -A d x +5 — 164 d x -2 d x +z + 356 d^-i d x +2 + 228 rf^+i) d y +i d y +A 
+ (—44 da;-4 d x +5 — 52 oii+3 + 228 d x -i d x +2 + 324 d x d x +i) d y +2 d y +z 

appears to vanish. 

Messy doubling formulae, in case I don't find anything better in the next few days: 



a, [ x + i> 

V- 3, 



d—y+ X +\ dy-\-a 
d-y-1 dy-2 
d\ — y dy 
di — y dy + 3 
d5-y dy + 4 



d2k-i = —{dk-5 dt-2 dk+A dk+7 + 3 dk-3 dk-2 dk+2 dk+7 — 7 d 2 k _ 2 d k +\ dk+7 
— 11 rffc-2 dk-i dk dk+7 + 3 dk-5 dk dk+A dk+5 + 9 dk-z dk dk+2 dk+5 — 41 dk-2 dk dk+i dk+s — 13 rffc-i d\ dk+5 
— 7 dk-s dk+i dl+A — 11 dk-5 dk+2 dk+z dk+A — 41 dk-z dk+i dk+2 dk+A + 89 dk-2 d\ +1 dk+A + 57 dk-i dk dk+i dk+A 
— 13 dk-z d\+2 dk+z + 57 d k -2 d k +± d k +2 dk+z + 81 d k -i d k d k +2 dk+z) /40 



d2k — (dk-6 dk-i dk+5 dk+s — 3 dk-A dk-i dk+z dk+s + H dk-2 dk-i dk+i dk+s 
-47dl_i dk dk+s + 3 d k -e d k +i dk+5 dk+e — 9 dk-A dk+i dk+z dk+e + 53 dk-2 d't+i dk+e 
-161 dk-i dk dk+i dk+e — 7 dk-e dk+2 d\+ 5 — 11 dk-e dk+z dk+A dk+5 + 41 dk-A dk+2 dk+z dk+5 
-137 dk-2 dk+i dk+2 dk+5 + 329 dk-i dk dk+2 dk+5 + 13 dk-A d\+z dk+A — 81 dk-2 dfe+i dk+z dk+A 

+577 d k -! d k dk+z d fc+4 )/120 



d2 k = {dk-e dk-i dk+5 dk+s — 3 rffc-3 dfc+2 dk+s + 8 dk-2 dk-i dk+i dk+s 
—50 dk-idk dk+s + 3 dk-e dk+i dk+5 dk+e ~ 9 dk-z dk+i dk+2 dk+e + 44 dk-2 d\ +1 dk+e 
— 170 dk-idk d k +i d k +e — 7d k -e dk+2 d\ +5 — 11 d k -e dk+z dk+A dk+5 + 41 d k -z d\ +2 d k +5 
—96 d k -2 dk+i d k +2 dk+5 + 370 d k -\ d k d k +2 dk+5 + 13 d k -z dk+2 dk+z dk+A — 68 d k -2 dk+i d k +z dk+A 

+590 dk-i d k dk+z dfc+4)/120 

As an existence proof, here is a 6 th order relation: 

= ( (dk-e dk-5 dk-z + dk-e d 2 k - b dk-4) dfc-i + (dk-e d\- h d k -z + dk-5 dk-A) d k -2 + d\_ 6 dl_ 3 

2 2 2 \ 2 / / 2 

+2d k -e dk-5 dk-A d k _ 3 + d k _ 5 d k _ A dk-z) d k + ((d fc _ 6 dk-5 dk-2 

2 2 \ 2 / 22/22 

+dk-e dk-A d k -z + 2 dk-e dk-5 d k -A ) d k -i + (2 dk-e d k -5 d k -2 + (2 d k -e dk-z 



26 



—56 4-6 4-5 4-4 dk-3 + 2 d k _ 5 4-4) dk-2 + 2 4-6 4-4 4-3 

+2 4-5 4-4 rffc-3) C?fc — 1 + 4-5 4-2 + (2 4-6 4-S 4-3 + 3 4-5 4-4 dfc-3) 4-2 

+ (2 4-6 4-4 4-3 + 2 4-5 4-4 dfe_a) 4 + (4-6 4-4 4-2 + 4-6 dfc-4) 4-1 

+ ((4-6 4-3 + 2 4-6 4-5 dk-ij 4-2 + (3 4-6 4-4 dk-3 + dk-5 4-4J 4-2 

+4-4 dk-3) 4-1 + ((24-6 4-5 4-3 + dL,5 4-4) 4-2+ ( 2 4-6 4-4 4-3 
+3 4-5 4-4 4-3) 4-2+2 4-4 4-3 4-2) 4-1 + 4-5 4-3 4-2 
+2 4-5 4-4 4-3 4-2 + 4-4 4-3 4-2- 

It appears that 

S—i(x), sq(x), = — 1,0, 1,1,1,1,1,— 2,x,x — l,2x — 3,x — 4, x 2 — 4x + 2, x 2 — x — 1,..., 

gives a (non-divisibility) sequence of polynomials. The special case x = 1 appears to give eight interlaced arith- 
metic^) sequences: 

s 89+ ,(l) = (-)«[0,l,2g + l,l,l,l,-2g-2,l] <r<7. 

A curious initialization is 

. . . , — ^=,1,1,1,1,1,1, -V2, 2 - V2, 3 - 2 V2, 5 - 4 V2, 10 - 8 V2, 28 - 20 V2, 107 - 76 \/2, 455 - 322 V2, . . . , 
v2 

where the first term with denominator > 1 is 

, 510156039514521981558192050 - 360734795003990787362927953 y/2 
d34 - , 

and the first term (after d-i — — 3/\/2) with magnitude > 2 is 

d 3 g « 2.3813134529. 

References: 

Plouffe's Inverter: http://www.lacim.uqam.ca/pi/indexf.html 

Sloane's Sequence Server: http://www.research.att.com/~njas/sequences/ 

Theta Functions in Macsyma [5] 

Zagier notes on Somos5 and elliptic curve: 

http://www-groups.dcs.st-andrews.ac.uk/~john/Zagier/Problems.html 
and g]. 

i91ike double sum involving seven empirical constants for Somos6: 
http://grail.cba.csuohio.edu/~somos/somos6.html . 

Jim Propp's Somos sequence site: http://www.math.wisc.edu/~propp/somos.html . 

General: Google search Somos sequence* and Elliptic divisibility. 
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3 Modular Theta Functions 



This section explores mod P analogs of some classical special functions. We're interested in the general problem of 
developing modular analogs for classical special functions of analysis. The modular versions of exponentiation and 
logarithms have been known for two centuries. These are easily generalized to modular trigonometric and hyper- 
bolic functions, and thence to elliptic functions such as sn, cn, dn, and Weierstrass p. The modular Gudermannian 
function gd converts between sin and tanh. A modular version of the amplitude function am converts between 
sin and sn. I've already mentioned modular polylogarithms. This note introduces modular theta functions. 



3.1 Elliptic Function Review 

The classical elliptic functions arose from trying to determine the arc-length of an ellipse, by integrating expres- 
sions involving the square roots of cubic or quartic polynomials. Elliptic functions are analytic, complex valued 
functions, that take complex number arguments. They have two periods. One period is usually taken to be a 
real number, and the other is necessarily complex. Elliptic functions have poles as their only singularities (at 
finite locations). The two periods make a period parallelogram, covering the complex plane in a regular pattern. 
The function values repeat in each parallelogram. There are two popular ways of discussing elliptic functions, 
depending on whether the basic integral is the square root of a cubic or quartic polynomial. The two ways are 
equivalent, but the choices lead to differences in details of the formulas. The Jacobian elliptic functions ([HE], 
chapter 16) lead more naturally to theta functions. The basic Jacobian elliptic functions are sn(«, fc), cn(u, fc), 
and dn(u,fc). Typically k is fixed in an application (it is related to the shape of the period parallelogram), and 
it is often elided to simplify formulas. These functions have two poles and two zeros in each copy of the period 
parallelogram. Some of the fundamental formulas are 

sn 2 u + cn 2 u = 1 
k 2 sn 2 u + dn 2 u = 1 



sn(u + v 



sn 2u 



sn u = cn u dn u 

sn u cn v dn v -(- sn v cn u dn u 



1 — k 2 sn 2 u sn 2 v 
2 sn u cn u dn u 



1 — k 2 sn 4 « 
sn u = sin(am u) 



The sn and cn functions can be regarded as sine and cosine of a distorted input. The am function captures the 
distortion. Chapter 16 of Abramowitz & Stegun 2, 3j has much more, and the classic Whittaker & Watson |llj 
explains what's going on. 

3.2 Elliptic Functions vs Elliptic Curves 

An elliptic curve can be parameterized by elliptic functions, just as a circle can be parameterized by circular 
(trigonometric) functions (cos and sin) and a hyperbola by hyperbolic functions (cosh and sinh). 

Y 2 = k 2 X i -{l + k 2 )X 2 + 1 
(X,Y) = (sn it, cnudnu) 
(cnudnti) 2 = (1 — sn 2 it)(l — fc 2 sn 2 w) 
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The original applications of these ideas were for real and complex applications, but there has been a minor 
component of number theory ever since Diophantus (c. 200 AD) introduced problems that reduced to cubic 
curves. The theory of finding integer and rational points on elliptic curves has undergone a major development 
in the last century. The application of elliptic curves to cryptography depends on the fact that what works for 
the fields R, C, and Q can often be made to work for the fields mod p and the Galois Fields GF[p fc ]. 



3.3 What's a Theta Function? 

Theta functions were probably first introduced by Euler. They arise from some infinite products related to the 
partition function. Theta functions have rapidly convergent series, and they are closely related to elliptic functions, 
which makes them useful in computing elliptic function values. 

Like elliptic functions, theta functions are complex valued functions with one complex argument, and a second 
shape parameter. They only have one true period, but they have a second quasiperiod. The period is often taken 
to be 1 or 2ir, while the quasiperiod is a complex number. Together the two define a parallelogram, as with elliptic 
functions. Changing the argument by one period leaves the value of the theta function unchanged, while changing 
the argument by the quasiperiod multiplies the theta function by a fixed value. In one sense, theta functions are 
easier than elliptic functions, since they have only one pole and one zero per parallelogram. There are 4 basic 
theta functions, with the pole located either in the corner of the parallelogram, in the middle of a side, or in the 
center. Elliptic functions are ratios of theta functions, with the same parallelogram. 

3.4 Some Basic Properties of Theta Functions 

oo 

V4{z,q) = } j (-) g e 



n— — oo 



4 (z + 7t) = 4 0) 

4 (z + T ) = e 2lT tf 4 ( 

01,2,3 = 04 + half periods tt/2, r/2, (tt + t)/2 

OO 

04(^,9) = ]\{l-q 2n ){l-q 2n - 1 cos2z + q in - 2 ) 



h , 03 
— , dn ~ — 



04 ' 04 ' 04 

Similar products exist for 0i,2,3- The partition generating function is l/Y[ =i(l — 9™)- Gosper has developed a 
package for manipulating theta functions in the symbolic algebra program Macsyma [5]- 



3.5 Theta Function Identities 

As with the elliptic functions, there is a multitude, nay, a plethora, of theta function identities. A small sampling 
is presented below. 0;(O) is abbreviated to 0j. 

0! + 0l = 0^ 

tflO)0l = 0304(z)-0l0 3 (2) 

Mv + *)Mv-*)<% = 0i(j/)04(2)-04(y)0i» 

4 (2y)0l = -*}(») 

0l(2y)0 2 0304 = 20! (2/)0 2 0)03 0)04(2/) 
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The equation #f + = 0| has no non-trivial rational solution. (Clearing denominators would lead to a solution of 
Fermat's equation for exponent 4.) This is in contrast to the situation with elliptic curves, where modular solutions 
of the elliptic curve equation ( "points on the curve" ) can be turned into rational solutions of an equivalent elliptic 
curve. 

The restriction against rational solutions can be circumvented by switching to another field where the equation 
has solutions. Two obvious choices are an algebraic number field, or to work modulo a prime number. I decided 
to experiment with the modulus P — 43. A 4k + 3 prime was selected so that more fourth powers exist and so 
that the number i = \J— 1 has independent meaning. I selected trial values mod 43 for 0i-4 and t?i_4(l). The 
values were chosen to be compatible with the equations above. Having selected these 8 numbers, we use theta 
function identities to compute a table of 0i_4 evaluated at 2, 3, etc. 

Theta Functions Mod 43 



N 


01 


02 


03 


04 








2 


14 


1 


1 


1 


11 


7 


2 


2 


11 


24 


33 


15 


3 


1 


1 


24 


24 


4 


34 


12 


32 


36 


5 


17 


25 


7 


3 


6 


35 





39 


30 


7 


38 


20 


3 


32 


8 


8 


25 


5 


11 


9 


16 


27 


40 


40 


10 


14 


32 


42 


23 


11 


24 


37 


39 


5 


12 





4 


15 


41 


13 


7 


34 


37 


29 


o 13/1 


7 


7 


-7 


-7 



Note that row 12 has a relationship to row 0, with i?i and 02 being doubled, while 03 and 04 are multiplied by -2. 
Rows 13 and 1 are similarly related, with ratios ±7. We find a ratio relationship with consistent signs, between 
rows 24 and 0, 25 and 1, etc.: 

Row 24 = 16 * row 0; row 25 = 24 * row 1; row 26 = 36 * row 2. 

With a large enough separation between rows, a constant ratio will develop, and after some larger separation, a 
true period will appear. 

An equation similar to the Somos recurrence works: 

0i(n) = 80i (n + l)0i (n - 1) + 0i(n + 2)0i(n - 2). 

We explore the ratio relationship between theta functions and elliptic curves: 

, v 03 0i (n) 
02 04 (n) 
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n sn(n) 

0-11 0,25,8,20,9,11,1,11,9,20,8,25, this row is symmetric 

12-23 0,18,35,23,... this row is the negative of the row above 

Cn and dn work as well. Checking out the elliptic curve that is parameterized by these modular elliptic functions, 

X — sn(n); Y — cn(n)dn(n) 

Y 2 = 1 - 7X 2 + 6X 4 = (1 - X 2 )(l - 6X 2 ) ; elliptic curve equation 
Cn and sn are relabeled sine and cosine: The set {cn(n) + isn(n)} is the powers of (13 — 2i) mod 43, but in an 
apparently random order. Note that 13 — 2i is on the mod 43 unit circle, since | j 13 — 2i j | = 13 2 + 2 2 = lmod 43. 
The random ordering is captured by am, which in the real world is the distortion function: 

sn(n) = sin(am(n)) 

3.6 Some Possible Uses of Modular Theta Functions 

Modular theta functions may be directly useful in Dime Hellman key exchange. They could be used to compute 
elliptic curve values. And they provide another example of a modular analog for a classical special function. 

4 Conclusion 

There's a lot more to learn here. 
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